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ABSTRACT 



Aims. We calculate bounds on the variation of the fine structure constant at the time of primordial nucleosynthesis and at the 
time of neutral hydrogen formation. We use these bounds and other bounds from the late universe to test Bekenstein model. 
Methods. We modify the Kawano code, CAMB and CosmoMC in order to include the possible variation of the fine structure 
constant. We use observational primordial abundances of D, ^He and ^Li, recent data from the Cosmic Microwave Background 
and the 2dFGRS power spectrum, to obtain bounds on the variation of a. We calculate a piecewise solution to the scalar field 
equation of Bekenstein model in two different regimes; i) matter and radiation, ii) matter and cosmological constant. We match 
both solutions with appropriate boundary conditions. We perform a statistical analysis using the bounds obtained from the 
early universe and other bounds from the late universe to constrain the free parameters of the model. 

Results. Results are consistent with no variation of a in the early universe. Limits on a are inconsistent with the scale length of 
the theory I being larger than Planck scale. 

Conclusions. In order to fit all observational and experimental data, the assumption I > Lp implied in Bekenstein's model has 
to be relaxed. 

Key words, early universe - cosmology: theory - cosmic microwave background 



1. Introduction 

Time variation of fundamental constants has been an ac- 
tive field of theoretical and experimental research since 
the la r ge nu mber hypothesis (LNH) was proposed by 



19831: iGleiser fc Tavloi^ ll985l : lOverduin fc WessonI Il997l) . 
where the gauge coupling constants may vary over cosmo- 
logical time scales. 



Dirad (|1937f ). The great predictive power of the LNH 
induced a large number of research papers and sug- 
gested new sources of variation. Among them, the at- 
tempt to unify all fundamental interactions resulted 



Following a different path of research, iBekenstein 



string derivec 


field theories ( 


Wu & Wand 19861 Maeda 


1988: 'Barr & MohaDatralll988; Damour & Polvakovll994l; 


Damour et al. 


,2002allbr). related brane- world theories 


(YoumI 2001allbl: 


Palma et al. 20031 


Brax et al. 


2003). and 


Kaluza-Klein theories (Kaluzalll921 


: Kleinlll926 


: Weinberd 
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(|1982[ ) proposed a theoretical framework to study the fine 
structure constant variability based on general assump- 
tions: covariance, gauge invariance, causality and time- 
reversal invariance of electromagnetism, as well as the 
idea that the Planck- Wheeler length (lO~'^^cm) is the 
shortest scale allowable in any theory, ft is well known, 
that bounds from the weak equivalence principle require 
I < Lp. However, in this paper we are going to analyze 
data from cosmological time-scales rather than planetary 
scales, the latter ones being relevant to probe the va- 
lidity of w eak equivalence prin ciple. The model was im- 
proved by I Barrow et al. diooi) using the main assump- 
tion that cold dark matter has magnetic fields dominating 
its electric fields. Moreover, a super s ymmetric general- 
ization of this model was performed by lOlive &: Pospelov 
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( 2002f ). allowing additional couplings between the scalar 
field and both a dark matter candidate and the cosmolog- 
ical constant. The model was also generalized in order to 
st udy the time var i ation of the strong coupling constant 
bv lChamoun et all (|200l[) . 

Different versions of the theories mentioned above pre- 
dict different time behavior of the gauge coupling con- 
stants. Thus, bounds obtained from astronomical and geo- 
physical data are an important tool to test the validity of 
these theories. In unifying schemes like the ones described 
above, the variation of each gauge coupling constant is re- 
lated to the others. In this paper, we limit ourselves to 
study the variation of the fine structure constant (a). 

The experimental research can be grouped in astro- 
nomical and local methods. The latter ones include geo- 
physical methods such as the natural nuclear reactor that 
operated about 1.8 x 10^ years ago in Oklo, Gabon, the 
analysis of natural long-lived P decayers in geological min- 
erals and meteorites and laboratory measurements such as 
comparisons of rates between clocks with different atomic 
number. The astronomical methods are based mainly in 
the analysis of spectra from high redshift quasar absorp- 
tion systems. Although most of the previously mentioned 
experimental data gave null results, evidence of time vari- 
ation of the fine structure constant was repo rted recently 
from high redshift quasar absorpt io n systems (IWebb et al 



obtain solutions for the early and late universe. In section 
IHwe show our results. Finally, in section [9] we discuss the 
results and summarize our conclusions. 

2. Bounds from BBN 

Big Bang Nucleosynthesis (BBN) is one of the most im- 
portant tools to study the early universe. The model is 
simple and has only one free parameter, the baryon to 
photon ratio rjB, which can be determined by comparison 
between theoretical calculations and observations of the 
abundances of light elements. On the other hand, data 
from the Cosmic Microwave Background ( CMB) provide 
an independent method for dete rmining r]B ( Spergel et al.l 
l2003i l200d ISanchez et alll2006f) . Considering the baryon 
density from WMAP results, the predicted abundances 
are highly consistent with the observed D but not with all 
■^He and ''Li. Such discrepancy is usually ascribed to non 
reported systematic errors in the observations of ''He and 
^Li. However, if the systematic errors of ''He and ^Li are 
correctly estimated, we may have i nsight into new p hysics 



beyond the minimal BBN model. iDmitriev et al.l (|2004l ) 



I999II2OOII: [Murphy et a"Lll2001bl l3. [200I iLevshakov et al 



20071 ). Ho wever, other recent independent analysis of sim- 
ilar data dMartmez Fiorenzano et al.l 2003: Q uast et al. 



Srianand et al.ii2004i) found no been eval uated by [B ergstro m et al 



2004 iBahcall et al.ll20ol ^ 
variation. 

Bounds on the variation of a in the early universe 
can be obtained using data from the Cosmic Microwave 
Background (CMB) radiation and from the abundances 
of light elements. Although these bounds are not as strin- 
gent as the mentioned above, they are important be- 
cause they refer to a different cosmological epoch. Finally, 
other bounds at redshift greater than 30 could be ob- 
tained from the 21 cm si gnal once it could be measured 
( Khatri fc Wandeltl l2007h . In this Daper, we perform a 
careful analysis of the time variation of a in the early uni- 
verse. First, we use available abundances of D, '*He and 
^Li and the latest data from the CMB to put bounds on 
the variation of a in the early universe without assum- 
ing any theoretical model. Then, we use these bounds and 
others from astronomical and geophysical data, to test 
Bekenstein theory. 

In section [2] we use the abundances of the light ele- 
ments to put bounds on where ao is the actual value 
of a, allowing the baryon to photon density r/B to vary. 
We also calculate the time variation of a keeping rjB fixed 
to the WMAP estimation. In section [3] we use the three 
year WMAP, other CMB experiments and the power spec- 
trum of the 2dFGRS to put bounds on the variation of 
a during recombination, allowing also other cosmological 
parameters to vary. In sections [U [5] and [6] we describe the 
astronomical and local data from the late universe. In sec- 
tion [7] we describe Bekenstein model for a variations and 



consider the variation of the deuterium binding energy in 
order to solve the discrepancy between D, '^He and ^Li 
abundances. 

In this section we focus on the possibility that the 
fine structure constant may be different from its ac- 
tual value during BBN. The dependence of the pri- 
mordial abundances on the fine struct ure co nstant has 

(I1999I) 



, 3^^^^^ . . and im- 
proved by iNollett fc Lo pez (200^ Semi-analytical anal- 
ysis have been perform ed by some of 



lier works (Landau et al. 200' 



oy some or us m ear - 
i IChamoun etldl l2007f ). 



Ichikawa fc Kawasakil ( 20021 ) studv the effects of varia- 



tion of fundamental constants on BBN in the context of 
a dilaton superstring model. In a following work, they 
study the primordial abundances of light elements when 
the fine structure constant and the cosmic expansion rat e 
take non-standard values ( Ichikawa fc Kawasakil 2004 ). 
Miiller et al. ( 2004 ) calculate the primordial abundances 



as a function of the Planck mass, fine structure constant, 
Higgs vacuum expectation value, electron mass, nucleon 
decay time, deuterium binding energy and neutron-proton 
mass difference and study the dependence of the last three 
quantiti es as func t ions o f the fundamental coupling and 
masses. Coc et al. 1 (I2OO7I) set constraints on the variation 
of the neutron lifetime and neutron-proton mass differ- 
ence using the primordial abundance of ^He. Then, they 
translate these constraints into bounds on the time varia- 
tion of the Yukawa cou plings and the fine structure con- 
stant. |Cxburt_et_alJ (j2005i) study the number of relativistic 
species at the time of BBN and the variation of fundamen- 
tal constants a and Gn and set bounds in these quanti- 
ties using the primordial abundances and the results of 
WMAP for rjB. 

In this work, we modify numerical code of Kawano 
()Kawanolll988lll992h in order to allow a to vary. In addi- 
tion to the dependences on a discussed by other authors. 



Mosquera et al.; Variation of a and the Bekenstein model 



3 



we als o include the depend ence of the light nuclei masses 
on a (jLandau et al.ll200a ). The code was also updated 



with the reactions rates reported by iBergstrom et al 

We consider available observational data on D 
^He and ''Li. For D we consider the values re- 

by 



i iorted 
200X1) ■ iKirkman et al 
1998aibl); ICrighton gtal 
[Oliveira et al . 
For ''Li ■? 



Pettini fc Bowen (120011) : 



O'Meara et al 



Buries &: Tvtler 



I 2004h : lO'Meara et al.l 1^200^ 



(2006). 
e consider 



200C); 'Bonifacio et al. 



t he re sults 
("19971) 



from I Ryan et al. 



Bpnifacio & Molard 
19971): ,Bonifac io ct al (2002D: lAsplund et al.. ()2006i) : 



Boesgaard et al.l (I2OO5I ): iBOTdfacio et al.l (|2007l) 



The He available observations can be summa- 



rized in the results reported bv iPeimbert et al 



Izotov et al.l (|2007h : lOlive fc SkillmanI (|2004h 



(120071); 
The re- 
ported values of ^He depend on t he ad opted set of Hel 



emissivities. In fact, llzotov et al.l (|2007|) report two val 



ues, o ne calculated with old atomic data (iBeniamin et al 



2002) and th e other with ne w atom ic data (| Porter et al 
20051 ) while iPeimbert et al.l (|2007( ) used the new val 



ues 
data. 



We consider the results ca lculated using new atomic 
Olive & Skillman (2004) reanalyze the values of 



Izotov fc Thuan (1998 ): Peimbert et al.. (,2000 ) for the pri- 
mordial abundance of ''He. They examine some sources 
of systematics uncertainties and conclude that the obser- 
vational determination of primordial helium abundance 
is limited by systema ti cs errors. We consid er the data of 
Peimbert et ail (|2007t ): 



Izotov et al.l (|2007f) in our analy 



sis. In table |T] we show the theoretical predictions of the 
abundances in the standard model (without variation of 
a) with rjB fixed to the WMAP estimate. 

Table 1. Theoretical abundances in the standard model 



Nucleus 


Our Code 


D 


2.569 X 10"^ 


*He 


0.248 




4.514 X 10-^° 



In order to check the consis tency of the data, we follow 
the analysis of Yao et al. ( 20061 ) for the data set considered 
in this work. We find that the ideogram method plots are 
not Gaussian like, suggesting the existence of unmodelled 
systematic errors. We take them into account by increas- 
ing the standard deviation by a factor S. The values of 
S are 2.10, 1.40 and 1.90 for D, ''He an d ^Li respectively. 
A scal ing of erros was also suggested bv lOlive fc Skillman 
(|2004 . 

We compute the light nuclei abundances and perform 
a statistical analysis to obtain the best fit values for the 
parameters, for two different cases: 

— variation of a allowing 77s to vary, 

— variation of a keeping r/B fixed. 



Even though the WMAP data are able to constraint 
the baryon density with great accuracy, there is still some 
degeneracy between the parameters involved in the statis- 
tical analysis. For this reason, we allow the joint variation 
of baryon density and the fine structure constant to obtain 
an independent estimation for rjB- 

2.1. Variation of a and rjB 

In this case, we compute the light nuclei abundances for 
different values of 775 and and perform a statistical 
analysis in order to obtain the best fit values for these pa- 
rameters. As pointed out by severa. 1 authors ( Cuoco et al. 
2004tlCvburtll2004l : [Co"c et al.ll2004l : llchikawa fc Kawasaki 
2002 . 2004h . there is no good fit for the whole data set 
even for — ^ 0. However, reasonable fits can be found 
excluding one group of data at each time (see table [2|) . 

For D-l-^'He, the value of 77^ is coincident with WMAP 
estimation and there is no variation of a within 3cr. On 
the other hand, the other groups of data, favour values 
far from WMAP estimation, and for D -I- ^Li, the result is 
consistent with variation of a within 6(t. In figures [T] to |3| 
the confidence contours and 1 dimensional Likelihoods are 
shown, considering the available data of ^He-|-^Li, D-|-^Li 
and D -|- ''He, respectively. 




3 3.1 3.2 3.3 3.4 3.5 3.6 
nB[io-'°l 



0.07 -0.06 -0.05 -0.04 -0.03 -0.02 



Fig. 1. Icr, 2(7 and Scr likelihood contours for ^ vs rjB 
(in units of 10"'") and 1 dimensional Likelihood using the 
data of ''He + ^Li. 



Table 2. Best fit parameter values and la errors for the 
BBN constraints on and jjB (in units of 10^'°). 



Data 


VB [10-"'] 




N-2 


D -h "He ''Li 


4.1881^^5 


-o.oosinS^ 


10.33 


*He + 'Li 


•J-^o»_0.148 


-0.051 ± 0.009 


1.00 


D-h'Li 


' •'^"^-0.491 


0.165 ±0.012 


1.00 


D-h^He 


R iqc,-f0-443 
D.iyO_Q 4]^g 


-0.019 ±0.008 


1.00 
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" 6 6.6 7 7.5 8 8.6 
ilB [10-'°1 



Fig. 2. la, 2(7 and 3cr likelihood contours for — vs 775 



(in units of 10 
data of D + ^Li 



10\ 



and 1 dimensional Likelihood using the 




0.04 -0.03 -0.02 -0.01 

Aa/an 



Aa 



Fig. 3. la, 2a and 3a likelihood contours for ^ vs 7]b 
(in units of 10^^") and 1 dimensional Likelihood using the 
data of D + ''He. 



2.2. Variation of a with i]b fixed 

Once again, we compute the light nuclei abundances for 
different va lues of k eeping 773 fixed to the WMAP 
estimation ( Spergel et al.ll2006) . and perform a statistical 
analysis in order to obtain the best fit value for — . As 

. . . . ao 

pointed out in the previous section, there is no good fit 
for the whole data set and for ''He + ''Li, even for — ^ 
(see table [31 . 

For D + ^Li, the result is consistent with variation of a 
within 6(7, meanwhile for D+'*He there is no variation of a 
within 3(7. In figures H] and [5] the 1 dimensional Likelihood 
is shown, considering the available data of D + ^Li and 
D + "^He, respectively. 

In order to test Bekenstein model, we consider the re- 
sults obtained using D + ''He in this section. 

It is worth to mention that if we don't consider ^Li data 
the results with and without varying rjB are the same. 



Table 3. Best fit parameter values and la errors for the 
BEN constraints on — . 



Data 


Aa 

an 


N-l 


D + "^He + ^Li 


0.077 ± 0.001 


23.28 


"He + ^Li 


0.077 ± 0.001 


45.18 


D-f ^Li 


0.129 ±0.006 


1.58 


D + ^He 


-0.020 ±0.007 


0.90 




0.1 0.12 0.14 0.16 0.18 0.2 

Aa/ao 

Fig. 4. 1 dimensional Likelihood of — using the data of 
D + ^Li. 




Fig. 5. 1 dimensional Likelihood of — using the data of 
D + ^He. 



3. Bounds from CMB 

Cosmic Microwave Background (CMB) radiation provides 
valuable information about the physical conditions of the 
universe just before decoupling of matter and radiation, 
and thanks to its dependence upon cosmological param- 
eters, it allows their estimation. Any change in the value 
of the fine structure constant affects the physics during 
recombination, mainly the redshift of this epoch, due to 
a shift in the energy levels and in particular, the binding 
energy of Hydrogen. The Thompson scattering cross sec- 
tion is also changed for all particles, being proportional 
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to o? . Therefore, the CMB power spectrum is modified 
by a change in the relative amplitudes of the Doppler 
peaks, and shifts in their positions. On the other hand, 
changes in the cosmological parameters produce simi- 
lar effects. Previous analysis of the CMB data includ- 
ing a possible varia t ion of a haye bee n performed by 
Martins et all (|2002h : iRocha et all lioosl ) : llchikawa et al ' 



In this paper, we use the WMAP 3-year tem- 
perature and temp erature-polarization power spectrum 



( Spergel et al.ll2006l). and ot her CMB ex periments such a s 
CBI (iReadhead et al.l|200l. ACBAR (IKuo et al.ll2004l) 



and BOOMERANG (jPiacentini et al.l 120061 : iJones et al 
2006 ), and the power spectrum of the 2dFGRS (jCole et al 



2005f ). We consider a spatially- flat cosmological model 
with adiabatic density fluctuations. The parameters of our 
model are: 



Aa 

ao 



(1) 



where ^CDAih^ is the dark matter density in units of the 
critical density, @ gives the ratio of the comoving sound 
horizon at decoupling to the angular diameter distance to 
the surface of last scattering, r is the reionization optical 
depth, Hs the scalar spectral index and As is the amplitude 
of the density fluctuations. 

We use a Markov Chain Monte Carlo method to ex- 
plore the parameter space because the exploration of a 
multidimensional parameter space with a grid of points 
is computationally pro hibitive. We use the p ublic avail- 
able Cos moMC code of Lewis fc Bridld ( 2002|) which uses 
CAMB (iLewis et al.ll2000l) and RECFAST (ISeaeer et al 



to compute the CMB power spectra, and we have 
modified them in order to include the possible variation 
of a at recombination. We ran ei ght different chains. We 
used the convergence criterion of Rafterv fc Lewi^ (1992) 
to stop the chains when R—l< 0.003 which is more strin- 
gent than the usually adopted value. Results are shown in 
table [4] and figure [6l Figure [6] shows a strong degeneracy 
between a and O, which is directly related to Hq, and 
also between a and Qsh^- The values obtained for ri^/i^, 
ft,, flcDAih^, T, and Ug are in agreement, within la, with 
the respective values obtained w ithout including any vari- 

Del). Our 



2ooe 



ation of a by the WMAP team (jSpergel et al. 
results are consistent within 2cr with no variation of a at 
recombination. 

We have also performed the analysis considering only 
CMB data. In that case, the strong degeneracy between 
a and Hq made the chains cover all the wide Hq prior, 
making it impossible to find reliable mean values and er- 
rors. Hence, we added a Gaussian prior to Hq, which was 
obtained from the Hu bble Space Telescope Key Project 
(jFreedman et al.ll2001r ). and chose the values of the mean 
and errors as those inferred from the closest objects in 
that paper, so we could neglect any possible difference be- 
tween the value of a at that redshift and the present value. 
In this way, we post-processed the chains, and found that 
the most stringent constraints were obtained in the first 
analysis (see figures [7] and [8|. 



Table 4. Mean values and errors for the principal and 
derived parameters including a variation. 



Parameter 


Mean value and Icr error 




n 021 ^+0 0009 




U.1UZ_Q QQg 


e 




r 






n ni 


71s 


965+° "" 


As 


3.039in^t 


Ho 


67.7ttl 



s 
< 




0.018 



0.02 0.022 0.024 



Fig. 7. 1(7 and 2a confidence levels contours obtained with 
CMB data with and without data of the 2dFGRS power 
spectrum. 




-0.08 -0.06 -0.04 -0.02 

Aa / an 



0.02 0.04 



Fig. 8. la and 2a confidence levels contours obtained with 
CMB data with and without data of the 2dFGRS power 
spectrum. 



4. Bounds from Quasar Absorption Systems 

Quasar absorption systems present ideal laboratories to 
search for any temporal variation in the fine structure 
constant. Quasar spectra of high redshift show the ab- 
sorption resonance lines of the alkaline ions like CIV, 
Mgll, Fell, SilV and others. The relative magnitude of 
the fine splitting of resonance lines of alkaline ions is 
prop or tional to a^. Seve r al aut h ors dCowie fc Songaila 
19951: IVarshalovich et all Il99d iMurphv et al.l l2001c : 
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Fig. 6. Marginalized posterior distributions obtained with CMB data, including the WMAP 3-year data release plus 
2dFGRS power spectrum. The diagonal shows the posterior distributions for individual parameters, the other panels 
shows the 2D contours for pairs of parameters, marginalizing over the others. 



Chand et aTlbOOsHMartfnez Fiorenzano et al.ll2003l ) stud- 
ied the SilV doublet absorption lines systems at differ- 
ent redshifts ( 2.5 < z < 3.33), to put bounds on the 
variat ion of a. iBahcall et al. (|2004|) used O III emission 
hues. IWebb et all (ll999^ : iMurphv et all (|2001bL ^qM) 
compared transitions of different species with far be- 
tween atomic masses and led to a single data consis- 
tent with time varying fine structure constant for a 
range of redshifts (0.5 < z < 3.5). Howev er, other re 



cent i ndependent analysis of similar data ([Quast et al 
2004 ISrianand et all 12004 iGrupe et all l2005f) found no 



variation. Another method, to test cosmologic al vari - 
ation of a, was proposed by iLevshakov et al. (|2005h 
from pairs of Fe II lines observed in individual expo- 
sures from a high resolution spectrograph. The authors 
found no variation of a at z = 1.15 and z = 1.839. 
However, a recent reanalysis of spectru m of the quasar 



QllO l-264 found variability within la (jLevshakov et al 



20071 ). We also consider in ou r analysis the boun d s men - 
tioned in I Wolfe et al.l (Il976l): ISpinrad fc McKed ([1979); 
Cowie fc Songaila (|1995^ : iTzanavaris et al.l (|2007l ). which 
obtained by con iparing the o p tical a nd radio 



were 



redshifts. Furthermore, iMurphv et al.l (|2001a[) compare 



molecular and radio lines and o btain a more stringent 
constraints. On the other hand, iDarling (|2004[ ) reports 



bounds on the variation of a at 2 = 0. 2467 from the 
satellite 18 cm OH conjugate lines. Finallv. [Kanekar et al 



(12005) compared the HI and OH main line absorption red- 
shifts of the different components in the z = 0.765 ab- 
sorber and the z = 0.685 lens toward B0218-I-357 to place 

/ 2 N 1-57 

stringent constraints on changes in F = ( ^ j 

Since we want to compare the prediction of a evolution 
with cosmological time, we consider each individual mea- 
surement of the papers cited above and not the average 
value reported in each paper. 

5. Bounds from Geophysical Data 

5.1. The Oklo Phenomenon 

One of the most stringent limits on time variation of the 
fine structure constant follows from the analysis of isotope 
ratios in the natural uranium fission reactor that operated 
1.8 X 10^ years ago at the present day site of the Oklo 
mine in Gabon, Africa. The proof of the past existence 
of a spontaneous chain reaction in the Oklo ore consists 
essentially of a substantial depletion of the uranium iso- 
topic ratio 235yy'238y -^i^ii respect to the current standard 
value in terrestrial samples and a correlated peculiar dis- 
tribution of some rare-earth isotopes. From an analysis of 
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nuclear and geochemical data, the operating conditions of 
the reactor could be reconstructed and the thermal neu- 
tron capture cross sections of several nuclear species be 
measured. 

The large values of the thermal capture cross sections 
of ^''^Sm,^^^Gd and ^^'^Gd are due to the existence of res- 
onances in the thermal region. In presence of such res- 
onance, the mono-energetic capture cross section is well 
described in the thermal region by the Breit-Wigner for- 
mula: 



t9 



{E- Er 



H 

4 



where p is the momentum of the neutron, g a statistical 
factor depending upon the spins of the compound nucleus 
of the incident neutron and target nucleus, Tj is the ra- 
diative width, r„ (E) is the neutron partial width, F is 
the total width. Thus, a shift in the lowest lying reso- 
nance level in ^^^Sm : A = ^i^^COkio) _ ^i49(now) 
derived from a shift i n the neutron capture cross section 

2OO0I : 



of the same nucleus (Fuiii et al 



Damour fc Dvson 



199Q). The shift in the re sonance energy can be translated 
( Damour fc Dvson 19961 ) into bounds on a possible differ- 
ence between the value of a during the Oklo phenomenon 
and its actual value. 

V arious authors (iDamour fc Dvson 1996t Fuiii et al 



2OOOI: iLamoreaux fc TorgersonI |2004 ) have analyzed the 



Oklo data in order to put bounds on a. iFuiii et al. 1 (|200d) 



used samples of ^^^Sm, ^^^Gd and ^^^Gd to reanalyze the 
bound on the resonance energy. They took the effect of 
contamination into account, ass uming the same contami- 
nation parameter for all samples. ILamoreaux fc Torgerson 
( 2004 ) employ a more realistic spectrum than the com- 
monly used Maxwell-Boltzmann to put the following 
bound on a: 



ao 



(-45± 15) X 10' 



(3) 



This bound is very similar to the one found by 
Fuiii et al. (l2000h and therefore we are going to consider 
it when testing the Bekenstein model. 

5.2. Long-Lived f3 Decayers 

The half-life of long-lived (3 decayers such ^^^Re has been 
used by several authors to find bounds on the variation of 
a. These nuclei have a very long half-life that has been de- 
termined either in laboratory measurements or by compar- 
ison with the age of meteorites. This last quantity can be 
measured from a decay radioactivity analysis. The most 
stringent bound on the variation of the half life, A, pro- 
ceeds from the comparison of ^^^ Re decay in the Solar 
System formation and the present (jOlive et al.ll2004r ): 



AA 



= (-0.016 ±0.016) 



(4) 



Sisterna fc VucetichI ( 1990l ) derived a relation between the 
shift in the half-life of long lived (3 decayers and a possible 



variation between the values of the fundamental constants 
(^i^QCD and Gi? at the age of the meteorites and their 
values now. They use a phenomenological model in which 
the abundance of any unstable nucleus will obey the fol- 
lowing decaying law: 



N ^ Nq exp 



\t + \t^/2 



(5) 



In this paper, we only consider a variation and therefore, 
the following equation holds: 



(2) 

A ao 



(6) 



where a = 21600 for ^^^Re. 



6. Bounds from Laboratory 

The comparison of different atomic transition frequencies 
over time can be used to determine the present value of the 
temporal derivative of a. Indeed, the most stringent limits 
on the variation of a are obtained using this method. The 
dependence of hyperfine transition frequencies with a can 
be expressed as: 



VHyp ~ — —RoocFRELiaZ) 
PB mp 



(7) 



where /i is the nuclear magnetic moment, fis is Bohr mag- 
neton, Roo is Rydberg's constant, uip and rrie are the pro- 
ton and electron mass and F^iel is the relativistic contri- 
bution to the energy. 

The comparison of rates between clocks based on hy- 
perfine transitions in alkali atoms with different atomic 
number Z can be used to set bounds on a*^^-^ where 
k depends on the frequencies measured and refers to 
the nuclear magnetic moment of each atom. The first three 
entries of table [6] show the bounds on ^ obtained by com- 
paring hyperfine transition frequencies in alkali atoms. 

On the other hand, an optical transition frequency has 
a different dependence on a: 



Vopt ^ RooBFi{a) 



(8) 



where B is a numerical constant assumed not to vary in 
time and Fi{a) is a dimensionless function of a that takes 
into account level shifts due to relativistic effects. The 
comparinson between an optical transition frequency and 
an hyperfine transition frequency can be used to set bound 
on a^^li^. 

nip fj,B 

Different 



auth o rs. | Bize et al. (Sooi); iFischer et al 



(|2004l ): iPeik et all (|2004j) . have measured different opti- 
cal transitions and set bounds on the variation of a us- 



ing different methods. IFischer et al.l (|2004f ) have consid- 



ered the joint variation of a and We have reanalyzed 



the data of IFischer et al.l (|2004l ). considering only a vari- 
ation, yielding the fifth entree of table jS) On the other 
hand, iPeik et al.l (|2004l ) have measured an optical transi- 
tion frequency in ^^^Yb with a cesium atomic clock. They 
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Table 5. The table shows the compared clocks, the value of ^ and its corresponding error in units of 10 °y r , 
the time inte r val fo r wh ich the variation wa s me asured and the referen ce. References (1) IPrestage et alj (Il99 5 |); (2) 



the time inte r val to r wh ich the variation wa s me asured ana the re tere nce. Keterences 1 1 1 [ rrestage et al.l IllUiJoll; ( 
Sortais et al] (|2000f) : (3) iMarion et all (|2003l ): (4) iBize et alJ (|2003f) : (5) iFischer et all (|2004l ): (6) iPeik et alJ (|20ol 



Frequencies 




At[yr] 


Reference 


Hg+ and H maser 


0.0 ±37.0 


0.38 


(1) 


Cs and Rb 


4.2 ±6.9 


2 


(2) 


Cs and Rb 


-0.04 ± 1.60 


5 


(3) 


Hg and Cs 


0.0 ± 1.2 


2 


(4) 


H and Cs 


1.14 ±2.25 


5 


(5) 


Yb and Cs 


-0.58 ± 2.1 


2.8 


(6) 



perform a linear regression analysis using this result to- 
gether with other opti c al tra n sition frequen c y mea sure- 
ments from iBize et al.l (|2003l ): iFischer et aP (|2004h . We 
have already considered the other data, therefore, we have 
reanalyzed the data, using only the comparison between 
Yb and Cs frequency, yielding the sixth entree in table [6] 

7. The Model 

In this section, we solve the equation of the scalar 
field, which drives the variation of a in the Bekenstein 
model. First, we obtain the analytical solution for 
Friedmann-Robertson- Walker (FRW) equation for two 
different regimes and assure continu ity of the solution 
and its derivative. Unlik e other works ( Barrow et al.l[20o5 
Olive fc Pospelovll2002l ) , we do not assume that the scalar 
field is connected with dark matter field. We consider the 
weak field approximation and so only the electrostatic con- 
tribution to the scalar field equation is relevant. In this 
framework, the electric charge can be expressed in the 
form: 



e = eoe(x'^) 



(9) 



being e a scalar field. The Lagrangian for a charged par- 
ticle of rest mass m and charge ege: 



where = L is Lorentz invariant and after a 

gauge transfor mation it changes only by a perfect deriva- 
tive. Following Bekenstein ( 1982t ). we obtain the Lagrange 
equations for ([T0|) : 



(10) 



d (mu^) eo 



(11) 



Sem + 



(13) 



where 



S, 



(14) 



^ IGttG . 

belongs to the gravitational sector of the theory, 

Sem - -(16^)-^ J F^^^ F^,{-gf'^d^x (15) 

is the electromagnetic action and 

Sr,^ = [ - ["™^^ + ieoeWA^] [x' - x'{t)] d-^xilG) 

is the matter action where the coupling between matter 
and the gauge field depends on the scalar field responsible 
for the variation of a. The action of the scalar field can 
be expressed as: 



5, = 



1 he 
2~P 



(17) 



where I is a scale length which is introduced due to dimen- 
sional reasons and one of the assumption of this theoretical 
framework \sl > Lp.\t will be shown (see apendix[X| that 
this latter condition implies violation of the weak equiv- 
alence principle. However, this requirement could be re- 
laxed, due to string theorie s considerations (iBachadbOOl 
lAntoniadis fc Piohnel[l999h . 

Varying the total action with respect to the gauge field, 
the modified Maxwell equations are obtained: 



where 



—u^^^S^ [x' - x' (r)] 



1 



And identify ^),^: 

F^, = e-^[{eA,),^-{eA^),,] (12) Olne 

Following iBekensteinI ( 1982 ). the total action can be where 
written as: 



C7 v-g 

as well as the equation of motion of the scalar field: 



he 



da 



[x' ~ x' (r)] 



(18) 
(19) 

(20) 
(21) 



Mosquera et al.; Variation of a and the Bekenstein model 



9 



In an expanding Un iverse and evaluati ng the r.h.s of equa- 
tion pp)) following iBekensteinI (jl982l ) , the following ex- 
pression can be obtained: 



dte 



nc 



(22) 



where ^ = is a dimensionless parameter which mea- 
sures the fraction of mass in the form of Coulomb en- 
ergy to the total m atter density pm- As suggested in 
Sandvik et alj (120021 ) we use ( = 10"'*. Integrating equa- 
tion ()22|) . in an expanding universe and using pm ~ 
we obtain: 



dte 
e 



8n U. 



ait) 



{t - Q 



(23) 



where tc is an integration constant and is the total 
matter density in units of the critical density (pc)- In or- 
der to solve the above equation we must first solve the 
Friedmann equation for the different regimes we are con- 
sidering. 

In a flat Friedman-Robertson- Walker (FRW) universe, 
the equation for the scale factor reads: 



dta 
a 



Oo 

ait) 



ait) 



-TIT ^ +f^A (24) 



with the initial condition a (0) — 0, and satisfying that 
a (io) = flQ = 1. In the above equation, we assume that 
the scalar field contribution is negligible. Usually, this con- 
tribution is proportional to (^) and we expect the vari- 
ation of a to be of order 10~^. 

The FRW equation has no analytical solution in terms 
of elementary functions when radiation, matter and cos- 
mological constant are considered. We build a piecewise 
approximate solution by joining solutions obtained by con- 
serving only some terms of the r.h.s of equation ([M)) . We 
solve the FRW equation for two different cases: a) radia- 
tion and matter and b) matter and cosmological constant. 
In such way, solution a) can be applied to nucleosynthesis 
and recombination of primordial hydrogen whereas solu- 
tion b) is proper for quasar absorption systems, geophys- 
ical data and atomic clocks. 

First, we integrate equation ([M)) considering only mat- 
ter and radiation. In order to get an analytical expres- 
sion for the scale factor as a function of time, we change 
the independent variable to confornial time r/ as follows: 
a_RA/d?7 — dt. Defining ^ = Harj, we can write: 



a-RM (C) 



The time can be expressed as follows: 



Horn 



12 



(25) 



(26) 



Now, we solve equation p4[) considering only matter 
and cosmological constant and obtain: 



aMcit) = 7^ 



n. 



1/3 



sinh ( -V^^A-ffo it - to) 



-arcsh 




-1 2/3 



(27) 



The expansion factor must be a continuous and smooth 
function of time, and in order to match both solutions, the 
following conditions have to be fulfilled: 

a_RA/(ii) — o-Mciti) 

^r^*^) = -^^^'^ 

from where, we obtain: 
aRAiiti) = aMciti) ' 



Ha 



1/4 



(28) 



In order to compare with astronomical and local 
bounds, we use the value of the cosmological parameters 
reported in Yao et al. ( 20061 ). 

Now we can solve equation (1231) using the equations 



3 and (EZI). Using In ^ ~ i^, we obtain the follow- 
ing expressions for the variation of a in the two different 
regimes. 

Defining A(^) = ^f^m + i^/Tl^ for t < ti: 

2 



Aa _ 1^ 

TT 



^1 In 



( m 



--cf- 



In 



1 



^mrr. (I 1 



A(a) 2vn: 



=ln 



m 
mi) 



aA(c) 



Ho (ti - tc) 



th ^arcsh 



_ Ho ito — tc) 

Defining v = Hgt, we can write for t > ti 

Aa 
ao 



3/4' 



(29) 



2tt 



V — Vr 



th {^^y/n^iv - uo) + arcshy'H^ 




^sinh QV^A iv - I'o) 



-|- arcsh 




(30) 



10 



Mosquera et al.; Variation of a and the Bekenstein model 



8. Results 

8.1. The early universe 

In section [2l we have used the primordial abundances of 
D, "^He and ''Li to put bounds on the variation of a in 
the early universe. First, we have performed a statistical 
analysis in order to check the consistency of each group 
of data and modified the observational errors accordingly. 
We have shown that all data could not be fitted at the 
same time, but reasonable fits can be found considering 
two groups of data at the time. We have analyzed the 
case where the baryon density is a free parameter and the 
case where it is fixed to the WMAP value. Tables and 
confidence contours are shown in section [51 In all analy- 
sis describe in section [2] (with or without allowing rjB to 
vary), we find that excluding the ^Li data, our results are 
consistent with WMAP estimation and no variation of a. 

In section |3] we have used the three year WMAP data 
together with other CMB experiments and the 2dFGRS 
power spectrum to put constraints on the variation of a 
during recombination. Tables and confidence contours are 
shown in section [3) 

We summarize our results of the variation of the fine 
structure constant in the table [S) Our results are consis- 
tent with no variation of the fine structure constant in the 
early Universe. 



Table 6. Best fit parameter values and la errors for the 
BEN and CMB constraints on 



Table 7. Best fit parameter values and la errors of the 
Bekenstein model. 



Group of Data 


Aa 
an 


BBN 


-0.020 ± 0.007 


CMB 


-0.015 ± 0.012 



8.2. The Bekenstein Model 

In this subsection we compare the Bekenstein model pre- 
dictions obtained in section [7] with astronomical and geo- 
physical data described in sections [H [5] and [H] and with the 
bounds on a from the early universe we have obtained in 
sections [2] and |31 

Fixing the time, equation ([59)1 or ([50)1 gives the predic- 
tion for a variation as a function of two free parameters: 

(^j-^ and Hotc. Therefore, we have N (number of 

data we are considering: 1 from Oklo, 1 from ^^''Re, 6 from 
atomic clocks, 1 from BBN, 1 from CMB, 274 from QSO) 
conditional equations with two unknowns. We perform a 

test to obtain the best values of the free parameters of 
Bekenstein's theory. Our results are shown in table [71 

We also have performed the same statistical analysis 
discarding bounds of each group of data. In most of the 
cases, the results are similar than those considering all 
data. However, discarding the bound from nucleosynthesis 



Data 






All data 


0.000 ± 0.003 


(3.2 ± 1.4) X 10"*' 


Without Oklo 


0.000 ± 0.014 


(3.2 ± 1.4) X 10"® 


Without ^**''Re 


0.000 ± 0.003 


(3.2 ± 1.4) X 10"^ 


Without atomic clocks 


0.000 ± 0.003 


(3.2 ± 1.4) X 10"® 


Without BBN 


0.000 ±0.017 


(3.4 ± 1.3) X 10"^ 


Without CMB 


0.000 ± 0.003 


(3.2 ± 1.4) X 10"® 



changes the value of i^Y^j Hotc several orders of magni- 
tude. Thus, the bound obtained from the primordial abun- 
dances of the light elements are crucial to fix the value of 

\ 2 

I 



Our results show that the available limits on a varia- 
tion are inconsistent with the scale length of the theory I 
being larger than Planck scale. 



9. Summary and Discussion 

In this paper, we have analyzed the variation of a in 
the early universe. We have modified the Kawano code, 
CAMB and CosmoMC in order to include the possible 
variation of a. We have used recent observational abun- 
dances of light elements to obtain bounds on — at the 
time of primordial nucleosynthesis. We have used recent 
data from the CMB and the 2dFGRS power spectrum to 
limit the variation of a at recombination. Results obtained 
in sections [2] and [3| are consistent with no variation of a 
during primordial nucleosynthesis and recombination of 
neutral hydrogen. 

It is important to check that the values of the baryon 
density obtained using the light elements abundances (sec- 
tion [5]) are consistent with the respective value obtained 
using data from the CMB (section [3]) . Using the relation 
rjB = 2.739x 10~*57s/i^, we find that results are consistent 
within la. 

We have also used our results from the early uni- 
verse and recent bounds from the late universe to test 
Bekenstein model. We have improved the analysis of 
the Bekenstein moclel wit h respect to a previous work 
( Landau &: VucetichI 12002 ) in various aspects: i) we have 
obtained analytical expressions for the Bekenstein model 
which include the dependence on tc (while other authors 
put tc = 0) for the variation of a in two regimes: a) radia- 
tion and matter and b) matter and cosmological constant, 
ii) the whole data set is updated. 
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On the other hand, E6tv6s-hke experiments provide 
stringent constraints on the Bekenstein model parameters. 
Constraints for j-^ can be set using these kind of exper- 
iments (see appendix [A|) . We obtain -j^ < 8.7 x 10"'^ 
which is one order of magnitude below the limits obtained 
in this paper using astronomical and geophysical data: 

< 6 X 10~^. Nevertheless, the importance of our anal- 
ysis lies on the fact that while Eotvos-like experiments 
test planetary scales, in this paper we test different time 
scales, namely cosmological time-scales. 

The values obtained for the free parameters of the 
model disagree with the supposition that I > Lp, im- 
plied in Bekenstein's framework. However, this latter re- 
quirement could be relaxed. Indeed, it should be noted 
that Bekenstein's framework is very similar to the dila- 
tonic sector of string theory, and it has been pointed 
out that in the context of string theories (|Bachas|[2000t 
Antoniadis &: Piolind 11999 1) there is no need for an uni- 
versal relation between the Planck and the string scale. 



Appendix A: Eotvos-like experiments 

A general expression for the Eotvos parameter and re- 
cent calculations on the prot on- neutron mass difference 
were performed by IChamoun fc Vucetich (2002.) . From 
this work it follows that: 



Aa 
a 



(A.l) 



where a is the acceleration for different bodies falling freely 
in a gravitational field g, Eem is the electromagnetic en- 
ergy of the falling bodies and M is the nucleon mass at 
rest. A bound on was estimated in the same paper: 
\Te\ < 1.2 X IQ-^. Comparin g this expression with equa- 



tion (45) of lBekensteinI (|1982D : 



a 2n \Lp ' 
it follows that: 



(A.2) 



(A.3) 



We obtain -p- 

Lp 



< 8.7 X IQ- 
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